Optical conductivity of Ortho-II YBa2Cu306.5 
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The Ortho-II phase of YBa2Cu306.5 is characterized by a periodic alternation of empty Cu and 
fiUed Cu-0 b-axis chains doubling the unit cell in the a direction. The extra oxygen in the full 
chains gives rise to an attractive potential for the holes in the planes. The planar bands split in 
two with a gap opening at the new Brillouin zone boundary k^, = ±7r/2, which we estimate from 
LDA calculations. Using a planar model which treats the d-wave superconductivity in a mean field 
approximation, we show that interband transitions produce a strongly anisotropic feature in the 
optical conductivity controlled by a region in fc-space close to (7r/2,7r/2). The edge position of this 
feature gives information on the temperature dependence of quasiparticle spectrum in this region. 
Bilayer splitting would show up as a double edge shape. 

PACS numbers: 74.72-h,74.25,78.30.-J,71.27+a 



Highly ordered Ortho-II YBa2Cu306.5, with an aver- 
age doping per planar Cu x ~ 0.1 and a Tc '^ 60 K, is 
a suitable system to study the properties of underdoped 
cuprates and the effect of the chains on the supercon- 
ducting planes in YBCO. The Ortho-II phase is charac- 
terized by a periodic alternation of filled and empty Cu-0 
b-axis chains, doubling the size of the unit cell in the a- 
direction"'^ . The oxygen ordering reduces the disorder and 
highly ordered samples have been prepared^. Unfortu- 
nately YBCO is not suitable for the angle-resolved pho- 
toemission (ARPES) measurements^ due to poor quality 
of the cleaved surface, but it is a convenient system for 
infrared experiments. 

In this letter we analyze how the oxygen ordering 
affects the band structure and optical conductivity of 
Ortho-II YBa2Cu3 06.5. The presence of oxygen in the 
chains gives rise to an attractive potential for the holes 
in the planes. Due to the alternation of empty and filled 
chains, this potential is 2a periodic. This causes a reduc- 
tion of the Brillouin zone in the GX direction (see no- 
tation in Fig. 1) with the corresponding splitting of the 
bands. We find a highly anisotropic optical conductiv- 
ity due to interband transitions between the split bands, 
characterized by a sharp-edge peak in the a direction, 
while it vanishes at threshold in the b direction. These 
features should be absent in non-Ortho II phases. The 
conductivity onset of absorption is controlled by a re- 
gion in fc-space close to {tt/2, 7r/2) and the position of 
the threshold depends on the quasiparticle dispersion in 
this region and on the value of the interband gap. 

The normal state of underdoped cuprates has an 
anisotropic gap, with an overall d-wave dependence. 
ARPES measurements suggest^ that with decreasing 
temperature, the opening of the gap removes larger por- 
tions of the Fermi surface (FS) resulting in Fermi arcs 
close to the nodal direction. The size of these arcs shrink 
to a point upon entering the superconducting state. In a 
bilayer material interplanar coherence appears as a split- 
ting of the planar bands, and consequently of the FS 
branches, into bonding (B) and antibonding (A) states 



eyi_s(k) = ep/(k) ± t±{'k). The existence of bilayer split- 
ting (BS) in the antinodal region in the overdoped regime 
is widely accepted^. On the contrary, results in the 
nodal directioni^, as well as in the underdoped and op- 
timally doped regimesiSiifliiLiSiiiii remain controver- 
sial. We show that the dependence of the Ortho-II inter- 
band transition threshold on the gap amplitude allows 
one to extract information on these questions from in- 
frared measurements. 

We have performed LDA^^ calculations using the TB- 
LMTO-ASA computation schemeiS. YBa2Cu306.5 crys- 
tallizes in orthorhombic structure (space grout Pmmm). 
The following lattice parameters obtained at room 
temperature^^ a=7.659A b=3.872l and c=11.725A were 
used in calculations The basis of muffin-tin orbitals con- 
tained (3s,2p,3d) orbitals for O, (4s,4p,3d) orbitals for 
Cu, (5s,5p,4d,4f) orbitals for Y and (6s,6p,5d,4f) orbitals 
for Ba atom. A k-point mesh of 6x4x2 was used for 
self consistency calculation loop. As suggested by re- 
cent NMR experiments^^, the 2a periodicity arising from 
the Ortho-II ordering produces a charge imbalance be- 
tween the Cu atoms below a filled and an empty chain. 
Evidence of charge modulation on the order of 0.05 in 
planar Cu has been found in resonant soft X-ray scatter- 
ing rateiS.. We have calculated the charge distribution in 
the Cu atoms using the same size for all Cu spheres. Our 
estimation for the charge modulation in Cu ions (includ- 
ing the occupation of all the orbitals in the basis) is on 
the order of 0.06 holes in the planes. 

The calculated band structure is shown in Fig. 1. The 
band marked with dots is the single chain band which 
crosses the Fermi level, associated with the filled chain. 
Rough estimate of the volume enclosed by the Fermi sur- 
face of the chain bands leads to 0.04 ± 0.01 holes per 
Cu in the planes. This value is smaller than expected, 
what could be due to deviations of LDA predictions close 
to the insulating phase. The effect of the 2a periodicity 
is reflected in the band structure. The planar B and A 
states each split in two bands. The splitting of the bands 
shows up clearly in the XS and RU directions, at which 




dispersion 



FIG. 1: Main figure: Band structure obtained from LDA cal- 
culations. Tlie band marlced witli dots is cliain-like. G — 
(0,0,0), X = (7r/2,0,0), S = (7r/2,^,0), Y = (0,7r,0), T = 
(0,7r,7r), R = (7r/2,7r,7r), U = (tt.O.tt), Z = (0, 0, 7r). Inset: 
LDA Fermi surface in tlie GXSY plane. The cliain brancli is 
enhanced with dots. Solid (dot-dashed) lines correspond to 
the split B (A) bands. 



the Brillouin zone has folded. The four planar bands 
appear in two pairs. Intra and inter-pair splitting are 
due to interlayer coherence and oxygen ordering respec- 
tively. Without the doubling of the unit cell both band 
pairs would be degenerate. The interband transitions 
discussed below connect lower (a) and upper (/3) band 
states with the same bilayer symmetry. LDA calcula- 
tions show a fc-dependent a{P) splitting. On the Fermi 
level, in the XS direction, it varies between 120 and 160 
meV independent of the bilayer band. The FS is shown 
as an inset. The chain branch is enhanced with dots. 
B and A branches are in solid and dashed-dot lines re- 
spectively. f3 bands have quasi-onedimensional character. 
This quasi-onedimensionality has been used to explain 
the anisotropy of the elastic scattering rate of Ortho-II 
YBa2Cu30ejj^°. B and A bands cut the XS line at Ba, 
Bj3 and Aa, Af3 respectively. Note that BS does not 
vanish at the nodes in LDA. 

To calculate the optical conductivity we use a planar 



e(k) = — 2i(cos k^ 



cos ky) + At' cos kx cos ky 

(1) 



2f"(cos 2kx + COS 2ky) — ^ 



kx,ky G (— 7r,7r) in lattice parameter units, which we 



take equal for a and b axes. The alternation of empty 
and filled chains is introduced by an external potential 
^Ejj,p=i,2,<Tc£j>cf,_^-,^ where i,j label sites in a and 
b direction, P the two superconducting planes. The in- 
terlayer hopping is tj^(k). In the reduced Brillouin zone 
kx e (-7r/2,V2) 
is 



ky G (— Ti", tt) the resulting dispersion 



eQ,/3(k) = - 2<cosA;j, - 2i"(cos2A;j; +cos2fcj^) -^±i^(k) 



± ( Acos'^ kx{t~2t' cos ky)'^ 



V^ 



2\ 1/2 



(2) 



First plus and minus signs correspond to L — A and B 
and second signs to a and /3, respectively. Since in LDA 
calculations V and tj_(k) are only slightly fc-dependent 
near the nodes, we set them constant. 

To describe the d-wave superconducting state we in- 
troduce a phenomelogical attractive potential of the form 
Vb(cos fca; — cos fcj,)(cos fc^ — cos ky) which reflects the sym- 
metry of the underlying Cu-0 plaquette. Thermal^i 
and microwave conductivity^^ measurements performed 
in Ortho-II YBa2Cu30e.5 have found that its super- 
conducting quasiparticles are well defined BCS excita- 
tions corresponding to an anisotropic low-energy gap 
with nodal lines. In the following, superconductivity is 
treated at the mean field level in the new split bands. 
Only intraband and interband pair scattering processes 
(7, k, a; 7, -k, a) -^ (7', k', a'; 7', -k', -a'), with 7, 7' = 
a, (3 are kept. The presence of the 2a periodic poten- 
tial breaks the rotation symmetry of the d-wave gap 
and, in the reduced Brillouin zone, it takes the form: 
AQ,_^(k) = iAo((7(k) cosfcj: — cosfcj^) with 



ff(k) = 



2 cos kx {t — 2t' cos ky ) 
'4cos2 kx{t-2t'co&kyY + ^f^ 



(3) 



The breaking of the tetragonal symmetry shifts slightly 
the position of the nodes, with a maximum shift at 
kx — 0. Here we take Aq = 71 meV for both B and 
A band3S.i2.2i. 

In linear response the real part of the optical conduc- 
tivity reads 



-"H = ?^E^^"W 



k,L 



e^(k)6^(k) + A„(k)A/3(k)\ 



EL(k)E^{]^) 






[1 - np {E^{k)) np (4(k))] 5 {u - (i?f (k) + ^(k))) 



e^(k)e^(k)+A^(k)A^(k) ' 
i?i(k)i?f(k) 



r 



[np {E^{\^)) - np (i?^(k))] 5 {lo - |i?^(k) - £;^(k)|) 



(4) 



Here, v and e are the volume of the sample and the 



A? «(k) are the 



electronic charge, E^^ = ^J <^'^}i{^) . ^„^ 

quasiparticle energies in the superconducting state, np 
is the Fermi distribution function and Af ** (k) the matrix 
elements in the a and b directions. 



M'"'(k) 



_ ¥"^{1-21' cos kyfsm^k^ 
4 cos^ kx {t 



M'"'(k) 






■ sm^ k„ cos^ fc^ 



4cos2fcx(t- 2t' cos kyY + "x 



(5) 
(6) 



The terms in parenthesis in |0J are the well known co- 
herence factors. When Ac((k) = A^(k) — the first 
(second) factor vanishes if the states involved in the tran- 
sition are on the same (different) side of the FS. The first 
term in Q accounts for the contribution of those tran- 
sitions which produce an excitation in both the a and 
(3 bands. The second term includes transitions which 
transfer a thermal excitation from one band to another. 

Upper part of Fig. 2 show the optical conductivity 
in the a and b directions for parameters which mimic 
reasonably well the LDA bands and FS in the GXSY 
plane. We take t = 360 meV, t' /t = 0.3, t" /t = 0.15, 
V = 120 meV ^i = -243 meV and t± = 135 meV. Fig. 
2a) show the a-axis conductivity in the superconducting 
(main figure) and normal, gapless, case (inset) at zero 
temperature. In the normal state the conductivity is 
characterized by a sharp single peak with an onset en- 
ergy w = 120 meV, while in the superconducting state 
it shows a double peak. At T = there are no excita- 
tions in the system and the second term in Q vanishes. 
Only those transitions which cross Ep are allowed. The 
minimum energy for the transition, and consequently the 
onset for absorption, is uj = V = 120 me\i^ and corre- 
sponds to the states with k^ — ±7r/2 between the two 
branches of the same bilayer FS, i.e. states between Aa 
and A(3, and between Ba and B(i (and the corresponding 
ones in the other quadrants) in the inset of Fig. 1 . Note 
the vicinity of these points to the nodal direction. Only 
these states contribute to the conductivity at threshold. 

A finite superconducting gap increases the energy of 
the transitions. The minimum energy still corresponds 
to a state at the edge of the zone, between the two nor- 
mal state FS branches. The double peak lineshape in 
the superconducting state marks the threshold energy for 
transitions in the B (lower onset) and A (higher onset) 
bands. The onset of interband transitions in the B-band 



occurs for fc-points close to the nodes in the supercon- 
ducting gap while in its A-band the relevant region of 
fc-space is further away from the nodes, leading to a dou- 
ble peak structure. The experimental observation of a 
double edge would be a clear evidence of the existence of 
BS in the underdoped regime. 

The peak lineshape is due to the k dependence of the 
matrix element M°°(k). M°'^(k) is finite and large at 
kx = i'''/2, and decreases away from the zone bound- 
ary, thus the conductivity decreases with increasing w. 
In the normal state, the peak in the conductivity is en- 
hanced by an approximate one-dimensional behavior of 
the transition energy for k between the two FS branches 
close to kx — ■n 12. This is also the reason of the larger 
magnitude of the first peak in the conductivity in the su- 
perconducting state. This strong enhancement would be 
smoothed by a ky dependent splitting at the XS direc- 
tion, as the one shown in Fig. 1, and, therefore, could be 
absent experimentally. 

Opposite to the behavior in the a-direction, M''''(k) 
vanishes at kx = ±7r/2, resulting in a vanishing a''^{uj) 
at the threshold frequency, as shown in Fig. 2b) for the 
normal (inset) and superconducting state. In the super- 
conducting state the BS splitting would appear in the 
b-axis conductivity as a double hump. The difference 
in intensity in the a- and b- directions is striking. The 
anisotropy in intensity and lineshape can be used exper- 
imentally to identify the interband transition. 

Correlation effects renormalize the band structure with 
respect to the LDA predictions. In the lower part of Fig. 
2 we plot the optical conductivity in the a direction, at 
zero (Fig. 2c) and finite (Fig. 2d) temperature, for the 
case in which the Fermi velocity has been reduced to 
the half of its previous value. The bilayer splitting has 
been adjusted such that the doping in the bonding (an- 
tibonding) band is 0.05 holes below (above) the average 
dopingii. Correlations are not expected to affect the 
Hartree potential due to the oxygen ordering, and we 
keep it unchanged. The qualitative features in the con- 
ductivity remain. Quantitatively the peak splitting has 
been slightly reduced and the bonding threshold deviates 
a bit from the gapless value. 

Fig. 2d shows the optical conductivity in the a direc- 
tion (7°° at T = 60K. At finite temperatures the existence 
of excitations reduce the contribution of the first term 
in (01) decreasing the amplitude of the zero temperature 
peak. The second term in|0J starts to contribute, and, in 
the presence of a gap, is responsible of a small structure 
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FIG. 2: Figs, a) and b) Optical conductivity in the a and 
b directions, respectively, in the normal gapless case (insets) 
and in the superconducting state (main figures) for the pa- 
rameters given in the text that mimic the LDA calculations. 
Note the different scale in the axis of a) and b). The axis and 
units of each inset are equal to the ones in the corresponding 
figure. Figures c) and d) show the the optical conductivity in 
the a direction in the superconducting state for the renormal- 
ized dispersion: t = 180 meV, /i = —156 meV and t± = 26 
meV, f'/f, t" /t and V are kept as a), see text. The values of 
/i and t± correspond to a doping xa = 0.15 and xb = 0.05 
holes in the A and B bands, respectively, c) shows the zero 
temperature conductivity for Ao = 71 meV. In d) the con- 
ductivity at 60 K corresponding to Ao = 71 meV (solid line) 
and to Ao = 35.5 meV (dashed line) is plotted. 



which appears at w ~ F and below. However, at the 
temperatures of interest T ^ Tc these effects are very 
small. On the other hand, if the superconducting gaps in 
the region of /c-space controlling these peaks are reduced 
with increasing temperature, the onsets for conduction 
and the peak splitting also decrease. This can be clearly 
seen in Fig. 2d, where the a-axis conductivities for two 
different values of Aq are compared. The evolution of the 
position of the peaks with increasing temperature gives 
us information on the modification of the superconduct- 
ing gap in the region of fc-space where the FS cuts the 
XS direction. 

In the pseudogap state, upon crossing Tc supercon- 
ducting coherence is lost, but a gap remains in the antin- 
odal spectrum. To describe the lose of coherence in (@J) 
we set AQ(k)A^(k) = in the coherence factors, but 
keep the same value (with gap) in the excitation energies 
_E^(k) and Ej^{k). This is equivalent to the condition, in 
the Nambu-Gorkov formalism, that the anomalous Green 
functions vanish. The energy of the interband transition 
does not change. Since in the superconducting state, the 
values of AQ,(k) and A^(k) are small in the area of fc- 
space which controls the conductivity close to threshold, 
this modification has little effect on the lineshape, and 
the discussion above remains valid even if superconduct- 
ing coherence is lost. 

In summary we predict an anisotropic feature in the 
infrared optical conductivity of Ortho-II YBa2Cu30(,,c,, 
which should be absent in non Ortho-II phases. It is 
a consequence of the oxygen ordering in every second 
chain. Possible bilayer splitting would cause a double 
edge shape. The threshold energy and its evolution with 
temperature give information on the quasiparticle spec- 
trum close to the nodes and on the effect of the chains 
on the planar electronic structure. Preliminary exper- 
imental results^'* have not shown any signature clearly 
identifiable with the one predicted here. The estimated 
intensity of this feature is comparable with the Drude 
contribution. Its appearance is based on the existence of 
quasiparticles. This picture seems to be well justified in 
the superconducting state at low temperatures and low 
energiesSiiSS.Its absence, if confirmed, could point out to 
a breakdown of the quasiparticle description at higher 
energies. 
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